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Abstract 

We use geometric measure theory to introduce the notion of asymp- 
totic cones associated with a singular subspace of a Riemannian manifold. 
This extends the classical notion of asymptotic directions usually defined 
on smooth submanifolds. We get a simple expression of these cones for 
polyhedra in E 3 , as well as convergence and approximation theorems. In 
particular, if a sequence of singular spaces tends to a smooth submanifold, 
the corresponding sequence of asymptotic cones tends to the asymptotic 
cone of the smooth one for a suitable distance function. Moreover, we ap- 
ply these results to approximate the asymptotic lines of a smooth surface 
when the surface is approximated by a triangulation. 

Introduction 

In the past decades, there has been a growing interest in defining geometric 
invariants describing singular spaces [2], pQ, [4], [5], [7], [8], [10], [TT], [13], [T7] . 
[19] . Such invariants are generally subject to two assignments: 

1. They must fit with the classical invariants when the underlying set is a 
smooth manifold or submanifold. 

2. They must satisfy some continuity conditions. For instance, if a sequence 
of singular spaces tends (for a suitable topology) to a smooth space, then 
the invariants defined on the singular spaces also tend to the smooth ones. 

After the length, the area and the volume, the most popular smooth geo- 
metric invariants are (sectional, Ricci, scalar, principal) curvatures, which are 
smooth functions (or tensors) defined on a (smooth) (sub) manifold and its tan- 
gent space. A classical approach to extending these curvatures to singular 
spaces W of a Riemannian manifold M, is to replace functions by measures 
on M. These measures are defined by integrating invariant differential forms 
over Borel subsets. As an example, Lipschitz-Killing curvature measures for 
singular spaces of a Riemannian manifold can be defined as follows : 
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• The first step is to generalize the unit normal bundle of a smooth subman- 
ifold to singular spaces. This has been done for convex subsets and for 
subsets of positive reach. More generally, the theory of the normal cycle 
[T8] allows us to define an integral current on a large class of singular sub- 
sets, called "geometric subsets" of a Riemannian manifold, generalizing 
the unit normal bundle. 

• The second step consists of defining standard differential forms on the 
tangent bundle of the ambient space. 

• Finally, integrating these differential forms on the normal cycle, builds 
invariant measures, satisfying the two assignments described above. 

We remark, however, that this building does not allow us to extract charac- 
teristic subsets induced by a couple (W, M) (where W is a singular space of a 
Riemannian manifold M) as asymptotic directions and principal directions. 

In the framework of geometric measure theory, more specifically using the 
normal cycle theory, we use the formalism introduced in [7], [8] and [17] to pro- 
pose an extension of the definition of the classical asymptotic directions to a 
large class of singular spaces W. Instead of building a new "curvature measure" 
on any couple (W, M), we associate a map that assigns a cone (called an asymp- 
totic cone) of vector fields on M to each Borel subset of M. In particular, if 
M = E N , we can reduce the target of this map to the set of quadratic cones of 
E^. Moreover, choosing a fixed scalar r > 0, we associate a field of cones leaving 
in the tangent bundle TK N with such a couple (W,'E N ). The two assignments 
described above are therefore satisfied. 

• The cones associated with a Borel subset reduced to a single point of a 
smooth surface lying in E 3 are nothing but the product of the classical 
asymptotic direction(s) by the normal line of the surface at this point. 

• If Wfc is a sequence of singular subsets (admitting asymptotic cones) , which 
tends to a smooth submanifold M, then the sequence of asymptotic cones 
associated with Wk built over a "regular" Borel subset B of M tends to 
the cone associated with M built over the same B. 

We give some applications of this building up for surfaces in E 3 . 

• If a plane field is defined on a singular surface of E 3 , we can define the 
asymptotic directions at a point m of the surface associated with a Borel 
set (for instance a ball centered at m) as the intersection of the asymptotic 
cone defined at m, with the plane defined at ra, to get two directions than 
can be called asymptotic directions. By integrating these directions, we 
get asymptotic lines depending obviously on the plane field. 

• Using a triangulated approximation of a smooth surface W 2 in E 3 , we 
build an approximation of its asymptotic directions at some fixed points 
by considering the intersection of the asymptotic cones of the triangulation 
with the planes spanned by the faces of the triangles. We can deduce 
an approximation of the asymptotic lines of W 2 if the triangulation is a 
"good" approximation of the smooth surface. 
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Finally, we mention that this approach can be extended to other character- 
istic subsets or subspaces, which could be the subject of future researches. 

The article is organized as follows: 

After a summary of the main notations pertaining to Riemannian submani- 
folds, the first section begins with smooth surfaces in E 3 . 

Since we define a new geometric invariant in quite a general context, our 
framework is an (oriented C°°) Riemannian manifold M (of finite dimension). 
After a review of the geometry of the tangent bundle of M (see [14] for in- 
stance), the second section introduces the so-called asymptotic form defined in 
the tangent bundle of M, easily derived from the fundamental form [8], [T7j . 

The third section uses the theory of the normal cycle [12], [18] to intro- 
duce the notion of asymptotic measures defined in the tangent bundle of M, 
associated with a large class of singular subsets Q (called geometric subsets in 
[12]). These measures are derived from the asymptotic form defined in Section 
[5] We describe them explicitly in classical situations (when Q is smooth or a 
polyhedron). Then, we construct the asymptotic cones. In particular, we show 
how this construction generalizes the asymptotic directions defined on a smooth 
surface in E 3 , fulfilling assignment [l] above. In particular, we give an explicit 
expression of the asymptotic cones associated with a 2-polyhedron in E 3 . 

The fourth section deals with a general theorem of convergence of asymptotic 
cones. This is the justification of our definition, and fulfills assignment [2] In 
particular, we give explicit assumptions in terms of the fatness of the sequence 
of polyhedra, implying the convergence of the sequence of cones for a suitable 
pseudo- metric. The last section presents two applications. First, we give a 
method to build asymptotic lines on a triangulation. Then, we give a method 
to approximate the asymptotic lines of a smooth surface approximated by a 
polyhedron. We test our method on various smooth or discrete surfaces. 

Some notation - For details on the theory of smooth Riemannian subman- 
ifolds, the reader may consult for instance, [6] or [9]. We use the following nota- 
tion. Let W n be an n-dimensional closed (oriented) smooth submanifold embed- 
ded in a smooth Af-dimensional (oriented) Riemannian manifold (M N , < ., . >). 
The manifold W n inherits a Riemannian structure by pulling back < ., . > by 
the embedding, which we still denote by < ., . >. We denote by TW n W n 
(resp. TM N ^> M N ) the tangent bundle of W n (resp. M N ) and E(W n ) 
(resp. E(M N )) the space of tangent vector fields over W n {resp. M N ). We de- 
note by V (resp. V) the Levi-Civita connection on (M N , < ., . >) (resp. W n ). 
We denote by T^W n ^? W n the normal bundle of the submanifold W n . The 
second fundamental form of the submanifold W n is the symmetric vector- valued 
(2, 0)-tensor 

h : TW n x TW n -+ T^W n 

defined as follows: 

Vx G TW n ^y e TW n ,h(x,y) = V x y - V x y. 
Let m be a point of W n . The isotropic cone 

CZ" = {xeT m W n :h m (x,x) = 0} (1) 
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of h m is classically called the asymptotic cone of W n at m. For any £ m G W n , 
the eigenvalues of < h m (.,.),£ m > are the principal curvatures of W n at m in 
the direction £ m . 

1 The case of smooth surfaces in e 3 

We first restrict our attention to smooth closed (oriented) surfaces W 2 embedded 
in the (oriented) Euclidean space (E 3 , < . >) bounding a domain D. Let 
£ be the normal vector field compatible with these orientations. The second 
fundamental form of W 2 can now be identified with the tensor < h (.,.)>£ > 
taking its values in C°°(W 2 ). We denote by Ai m ,A2 m the principal curvatures 
of W 2 at the point m, that is, the eigenvalues of /i m , by G the Gauss curvature 
of W 2 , that is, the determinant of h and by H its mean curvature, that is, its 
trace (in an orthonormal frame). In a frame of principal vectors (ei m ,e2 m ) at 
m (that is, eigenvectors of h m ) the matrix of h m is 



At each point m G with negative Gauss curvature, the asymptotic cone 
is the union of two lines. Integrating the corresponding vector fields gives 
rise to foliations of W 2 by the so-called asymptotic curves. By definition, at 
each point, the principal normal directions of these curves (considered as curves 
in E 3 ) are tangent to W 2 . 

Using measure theory, the goal of this paper is to define and study analogous 
cones associated with a (regular or) singular subspace of a Riemannian manifold 
M N , lying above any Borel subset of M N . 

Let us begin by explaining how we define such cones over any Borel subset of 
E 3 in the regular case; that is, when the subset is a (compact) domain bounded 
by a smooth surface W 2 . Let IV2E 3 be the tangent bundle of E 3 restricted 
to W 2 . If x is any vector field on 7V2E 3 , we build a signed measure ^>^ 2 as 
follows: For any Borel subset B of E 3 , we write 



where prT m w 2 denotes the orthogonal projection on T m W 2 and dm is the 



quadratic in x. If we force x to be a constant vector field, then we get a 
quadratic form (that we still denote by $* D {B)) on E 3 . This quadratic form has 
generically three eigenvalues, \\(B) , \2(B) , \$(B) , that we call the principal 
curvatures of B. The corresponding eigenvectors are called the principal vectors 
of 5, and the matrix of ^^(B) in this frame is 





(2) 




MB) 



0 



0 
0 



0 
0 



MB) 



0 



A 3 (B) 
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From this construction, we also introduce the isotropic cone associated with 
*'d(B): 

C D = {^E 3 ;^(B)=0}, (4) 

(the notation coming from the fact that we restrict our isotropic cone to constant 
- that is, parallel - vector fields in E 3 ). We call it the asymptotic cone of B (with 
respect to D). 

To clarify that this construction is linked with the classical pointwise situa- 
tion, suppose that B is reduced to a point {m} G W 2 . If y is a constant vector 
field such that y m e T m W 2 , then 

$ y D ({m}) h m (y m ,y m ). (5) 

If z is a constant vector field such that z m = £ m , then 

QUI™}) = 0. (6) 

This implies that in the frame (ei m , e2 m , £ m ), the matrix of ^>^({m}) is 

0 A 2m 0 . (7) 
0 0 0/ 

Consequently, the asymptotic cone C^ r,D is nothing but the cone spanned by 
the normal £ m and . 

This construction has some advantages: One can define the asymptotic 
cones at different scales by scaling the Borel sets (for instance, by taking balls 
of radius ^ as Borel subsets). Generically, we get three geometric invariants, 
Ai(.B), A2CB), A3 (£?), instead of two, and a two-dimensional cone instead of the 
union of two lines. Moreover, another important advantage of replacing func- 
tions by measures, is that this framework can be used for a large class of singular 
spaces (for instance, polyhedra, algebraic subsets, subanalytic subsets) of any 
codimension in any Riemannian manifold, as long as we can extend the notion of 
normal space. For instance, if one replaces the smooth surface by a polyhedron 
P bounding a domain D, we will get the follo wing explicit simple expression 



approximating the cone C^ r,T> (see 35 and also 33 for the exact formulas) 



$^(B)~^/(enB)Z(e)<i,e- > 2 , 

and 

C p B ar > V - {xGE 3 :^l(en B)Z(e) < x, e~ > 2 = 0}, 

(see Section [3] for the notation). This is why the theory of the normal cycle, 
extensively studied over the last decades [12], [18], will be our framework. In 
the following we describe the construction of asymptotic cones in such a large 
context. 
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2 Asymptotic forms 



To be self contained, we begin with a summary of the geometry of the tangent 
bundle of an (oriented) TV-dimensional smooth Riemannian manifold (M N , < 
.,. >). The reader may consult [T3] , [T5] and [8] for details. We denote by 

TTM N * T A N TM N the tangent bundle of the manifold TM N . As usual, we 
consider the exact sequence of vector bundles: 

0 — > TM N x MN TM N TTM N -A TM^ x m a, TM n — > 0, (8) 



where i denotes the natural injection defined by 

d 
dt 



i(u u u 2 ) ^(u! +tu 2 ) (9) 



and 

j = (K TM N,dit M N). (10) 

The vertical bundle of M is the subbundle V(M N ) = ker j of TTM N . The 
morphism i induces an isomorphism: 

i : TM N x MN TM N F(M N ). 

If m G and x G T m M N , the vertical lift of z e T m M N at x is the vector 
z v = i(x,z). The morphism J = i o j is an almost tangent structure on M N 
(J 2 = 0) and V(M N ) = ker J. Let 

5 : TM N — > TM N x m n TM N (11) 

be the canonical vector field defined by S(x) = (x,x) and let C : TM N — > 
V(M N ) be the vertical vector field associated with the (global) one-parameter 
group of homotheties with positive ratio, acting on the fibers of TM N . We have 
C = i o 5. We write n = pr 2 o i , where pr 2 denotes the projection on the 
second factor of TM N x m n TM N . Since M N is endowed with a Riemannian 
metric < •, • > and its Levi-Civita connection, we can build the corresponding 
right splitting 

7 : TM N x m n TM N — > TTM N (12) 

of the exact sequence j^J (satisfying j 07 = ld TM N x mnTM n). Let m G M N and 
x e T m M N . The horizontal lift of z e T m M N at x is the vector z h = j(x, z). We 
denote H X (M N ) = lm(j(x, •)), from which we construct the horizontal bundle 
H(M N ) such that, for all x G T m M N 

T X TM N = V X (M N ) © H X (M N ). (13) 

We denote by V(M N ) (resp. T~L{M N )) the space of vertical (resp. horizontal) 
vector fields. We denote by h : TTM N — > H(M N ) the horizontal projection, 
and by v : TTM N — > V(M N ) the vertical projection. We remark that h = joj. 
The morphism 

K = rj o v : TTM N TM N 

is the connector associated with the Levi-Civita connection. At every point 
x G T m M N , the morphism 

(d7T M x K) x : T X TM -+ T m M N x T m M N 
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is an isomorphism that identifies V X (M N ) with T m M N and H X (M N ) with 

T m M N . The bundle TTM N TM N is canonically endowed with the Sasaki 

metric < •, • > defined by the following conditions: 

' V(M N ) and H(M N ) are orthogonal, 

< i is an isometry, (14) 
7 is an isometry. 



If 



a:M N ^M N xM N 




is the diagonal map defined by a(m) = (m, m), then for every x G TM N , 
(d7TM x i^)" 1 o da(x) = x v © a/ 1 . 

Finally, the bundle TTM N TM N is also endowed with an almost complex 

structure F (F 2 = —Id), defined by the following conditions 

(15) 

Therefore F\ v : V(M N ) — > H(M N ) and F\ H : H(M N ) — > V(M N ) are 
isometries. In this Riemannian context, we give the following definition: 

Definition 1 1. The vector valued (N — l)-form on TM N defined for each 
X e H(M N ) by 

h x = [*Hod g e(i ? C A X)] A FX (16) 
is called the asymptotic (N — l)-form on TM N . 

2. The vector valued (N — l)-form on M N defined for each x G E(M N ) by 

h x = [* Uodge (FCAx h )]Ax v (17) 

is called the asymptotic (TV — l)-form on M N . 

In this definition, *Hodge denotes the Hodge duality on H X (M N ) for each 
x G TM N . (The introduction of the Hodge operator in the definition of the 
generalized second fundamental form can be found in [3 when the ambient 
space is Euclidean. We adapt it here to the general Riemannian situation. It is 
equivalent to the initial definition given in |5]). Using the identification of vector 
fields and 1-forms induced by the Riemannian structure, *Hodge(^C A X) is a 
(TV — 2)-form on H(M N ). On the other hand, FX is (identified with) a 1-form, 
null on H(M N ) and acting on V(M N ), and *Hodge (FC AX) is (identified with) 
an (N - 2)-form null on V(M N ) and acting on H(M N ). 



3 Normal cycles, asymptotic measures, asymp- 
totic cones 

3.1 Currents and normal cycles of singular spaces 

Let Vi(TM N ) be the space of l-currents of TM N (0 < I < 2iV); that is, the 
topological dual of the space V l (TM N ) of I -differential forms with compact 
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support on TM N , endowed with the topology of uniform convergence on any 
compact subset, of all partial derivatives of any order. The duality bracket will 
be still denoted by < •, • > if no confusion is possible. The space T>i(TM N ) is 
naturally endowed with the weak topology: if (Cfc)fceN * s a se Q uence of /-currents 
of TM N and if C is a /-current of TM N , then 

lim C k = C <=> Vcj G V\TM N ), lim < C^cj >=< C> > . (18) 

An /-current is rectifiable if it is associated with a rectifiable subset (see [17] for 
details). An /-current is integral if it is rectifiable and its boundary is rectifiable. 

When it exists, the normal cycle of a (compact singular) subset W of a 
Riemannian manifold M N is a closed integral current N(W) G T>n-i(TM n ), 
which is Legendrian for the symplectic structure on TM N dual to the canonical 
one on T*M N in the duality defined by the metric. The normal cycle is the 
direct generalization of the unit normal bundle of a smooth submanifold. Its 
formal definition was given in [12] . Although the normal cycle cannot be defined 
on any compact subset of M N , it exists for a large class of subsets, as convex 
subsets, polyhedra, subsets of positive reach, subanalytic subsets for instance. 
Following [12] , any compact subset Q of M N such that N(^) exists is said to 
be geometric, and N((/) is called its normal cycle. One of the main properties 
of the normal cycle for our purpose is its additivity [T2] : 

Proposition 1 // Q\ and Q 2 are geometric, then Q\ U Q 2 and Q\ n Q 2 are geo- 
metric and 

N(£i u g 2 ) = n(£i) + n(&) - N(& n g 2 ). (19) 

Here are some classical examples: 

1. The normal cycle of a smooth submanifold of a Riemannian manifold is 
the closed current associated with its unit normal bundle. 

2. If D is a compact domain whose boundary is a smooth hypersurface, then 
its normal cycle is the closed current associated with its outward unit 
normal vector field. 

3. If C is a convex body, then its normal cycle is the closed current associated 
with the oriented set 

{(m,0 :medC^e E 3 , ||£|| = 1,V* G C, < f , m2 >< 0}. 

4. The normal cycle of a polyhedron of K N can be computed by applying |l9| 
to a decomposition of the polyhedron into (convex) simplices and usingW 

3.2 Asymptotic measures, asymptotic cones 

Let us now define an asymptotic (signed) Radon measure on M N (resp. TM N ) 
associated with a geometric subset. We denote by B m n (resp. B tm n) the class 
of Borel subsets of M N (resp. TM N ) with compact closure. 

Definition 2 Let Q be a geometric subset of M N . 
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• The asymptotic measure defined on TM N , associated with Q and X G 
H(M) is the map 

®g ' B tm n — > R (20) 

defined as follows: 

\/B e Bt M n,*$(B) =< N(G), X Bh x > . (21) 

• The asymptotic measure defined on M N , associated with Q and x G 
S(M^) is the map <frg defined as follows: 

VB G B MN ,%{B) =< N(g), X „- 1{B) h* > . (22) 

If Q and B are fixed, the map 

x^<N(a),x OT -i (B) h*> 

is quadratic, inducing its isotropic cone. This remark leads to the following 
definition: 

Definition 3 Let Q be a geometric subset of M N . With any Borel subset B G 
B m n , we associate the cone 

C% = {xe E{M N ) : = 0} (23) 

and the cone 

Cg ar ' g = {x G E(M N ) : x parallel, = 0}, (24) 

which we call the asymptotic cones associated with B. 

In many applications, and for simplicity, it is easier to consider the cone of 
parallel vector fields, identifying a parallel vector field with its value at any point 
m G B. Obviously, C^ ar '^ C C£. We remark, however, that this new definition 
can be quite restrictive depending on the geometry of M N . For instance, if M N 
has non-zero constant sectional curvature, the only parallel vector field is the 
null vector field. In contrast, if M N = M N , the space of parallel vector fields 
is the space of constant vector fields, which can be identified with E^, which 
is much easier to manipulate. We deduce easily from [22] explicit expressions of 
these curvature measures in some particular cases. 

• The case of smooth submanifolds - Let W n be a (compact smooth) 
submanifold (with or without a boundary) embedded in M N and X G 
H(TM N ). Since the normal cycle of a smooth submanifold is its unit 
normal bundle, we deduce from [8] or [17] (Corollary 16 page 215) that 
for any B G B tm n, 

( B ) = / h^(pY TWn d7r M (X),pY TWn d7r M (X))d£dv, (25) 

Jsr ± w n nB 

where h£ denotes the second fundamental form of W n in the direction 
of the unit vector £, ST^W n denotes the unit normal bundle of W n 
and pr TWn denotes the orthogonal projection onto the tangent bundle 
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TW n . In particular, let W N ~ X be a (smooth oriented) hypersurface of 
M N bounding a domain D. (This assumption is not restrictive in our 
case, since our results are local. It allows to simplify some technical points 
by considering "only one side" of the normal cycles (the one corresponding 
to the outward unit normals), we have, for any B in B m n, 

®d( b ) = / h(pT TW N-idnM(X),pT TW N-idn M (X))dv, (26) 

Jw N ~ 1 nB 

where £ is the outward (with respect to D) unit normal vector field of 
W N ~ 1 , and h is the second fundamental form of W N ~ X in the direction 
£. We have then a correct generalization of [2] Consequently, 

— If B is reduced to a point m, 

c m n = i x e T mM N : h m (pr TW nx,pr TW nx) = 0}. (27) 

We deduce that C^ n is the cone spanned by C^ n and T^W n ; that 
is, we have the direct generalization of the corresponding cone defined 
for surfaces in E 3 in Section [11 

— If W N ~ X is a (smooth-oriented) hypersurface of E N bounding a do- 
main D, then 

C% = {x eTM N : / h(pr TW N-ix,pr TW N-ix)dv = 0}, 

Jw N ~ 1 nB 

(28) 

and 

£par,£> _ |^ ^ jgAT . / h(pr TW N-i X, pr TW N-i x)dv = 0}. 

Jw N ~ 1 nB 

(29) 



• The case of polyhedra - We will extend [28] and [29] to polyhedra. Let 
V be a domain in K N bounded by a (TV — l)-dimensional polyhedron 
pN-i p or ^ T-L(P N ), we can evaluate <£>^ above each simplex. In 
particular, if a N ~ 2 is a (TV — 2)-simplex, the support of N(V)\ a N-2 is the 
product Ccr, where C is a portion of circle. Let (ei, ...,ejv-2) be 

an orthonormal frame field tangent to a N ~ 2 . Any point of a N ~ 2 x C(j 
is a couple (m, ejv-i), where m is a point of a N ~ 2 and ejv-i is a unit 
vector orthogonal to a N ~ 2 . With these notations, we deduce from [8] or 
[T7] (Theorem 72 page 216) that for any 5 e B m n, 

E / <X,e^_ 1} > 2 . (30) 

We also deduce that for any 5 E Sj^jv, 

Cg*' 2 * = {xeE N : [ <x\ > 2 = 0}. (31) 

J(a N - 2 nB)xC 



31 



In particular, if V is a domain of E 3 bounded by a polyhedron P 2 , 
can be reduced to an explicit simple expression: First of all, we identify 
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the (vector) plane e x orthogonal to any (oriented) edge e of P with C, as 
follows : Let n\ G e x (resp. ri2 G e x ) be the unit (oriented) normal to 
the faces /i, (resp. $2) incident to e. Let e + (resp. e~) be the (oriented) 
normalized vectors spanning the bisectors of ri\ and ri2 (so that (e + , e~, e) 
is a direct frame of E 3 ). Any vector ae + + be~ of e 1 - is now identified with 
the complex number a + z6. An explicit integration over each term of type 



(a 1 HB) x C in 31 gives the following expression of the asymptotic measure 



and asymptotic cone: 



Proposition 2 1. For any B G B^s and any constant vector field x of 
= [(^(e)!^^!! 2 + sin(Z(c))7£((pr c ^a;) 2 )] , 

eGE 

(32) 

where E denotes the set of edges of P 2 , Z(e) the angle of the normal 
to the faces incident to e (being positive if and only if e is convex), 
pr e ± the orthogonal projection on e 1 - , and 7Z(pr e ±x) 2 the real part of 
the complex number (pr e ±x) 2 . 

2. In particular, 

= {x € E 3 : J2 l -^^[(Ae)\\pr e ,x\\ 2 + S HAe))n((pr e ,xf] = 0}. 

(33) 



We remark that if the Z(e)'s are "small enough" (this can happen for instance 
when P approximates "smoothly" a smooth surface), then sin(Z(e)) is "close 
to" Z(e) and 

®vi B ) ~^l(en B)Z(e) < x, e~ > 2 . (34) 

eGE 

And then, 

C p B ar ' V - {^E 3 :^l(en B)Z(e) < x, e~ > 2 = 0}. (35) 

eGE 

After choosing a scale r, we construct from the previous construction, a cone 
subbundle of TM N associated with a geometric subset of M N . Let us denote by 
B(m,r) the ball of radius r, centered at m G M N . With each point m G M N , 
and for a fixed (small enough) real number r > 0, we associate the cone C^f r y 

Definition 4 We call U tpnC^J? n the cone subbundle of TE N at scale r 

mG-Ui B(m,r) J 

associated with Q. 

We remark that the dimension of each fiber may change with m. This bundle 
is defined over the whole M N , even at the points m which are "far" from Q. 
If B(m,r) does not intersect the support of Q, then Cj^f r ) = T m M N . This 
phenomenon is visualized in Figure [l] where the cone la has its vertex (the black 



point) on the catenoid, and the cone [Tm has its vertex out of the catenoid. 
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(a) (b) (c) 

Figure 1: Some asymptotic cones built on a smooth surface (here on a portion 
of a catenoid) (a) An asymptotic cone in green whose vertex is at the center of 
its corresponding Borel set (a transparent ball) (b) The center of the ball may 
be out of the surface (c) When the Borel set is reduced to a single point on the 
surface, the asymptotic cone degenerates to the union of two planes 




Figure 2: Asymptotic cones built on a non-smooth surface, here a triangulation: 
the top of the Lilium tower (designed by architect Zaha Hadid) in Warsaw, 
Poland. 




Figure 3: Asymptotic cones built on a non smooth-surface. Here another tri- 
angulation: the Heydar Aliyev Center (designed by architect Zaha Hadid) in 
Baku, Azerbaijan 
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4 Convergence and approximation results 



In the previous paragraphs, we gave our first justification of the denomination 
of asymptotic cones. For surfaces W 2 in E 3 , this cone reduces at each point m 
of W 2 to (the product of the normal line by) the standard asymptotic directions 
of W 2 at m. We give now a second justification in terms of the convergence of 
sequences of polyhedra of E^. For simplicity, we restrict our study to (oriented) 
smooth hypersurfaces or polyhedra of K N bounding a (compact) domain. We 
will show, in particular, that, if a sequence of domains whose boundaries 
are polyhedra P/^ -1 , converges to a domain D whose boundary is a smooth 
hypersurface W N ~ 1 (in a sense that will be clarified later), then for a large 
class of Borel subsets P, the sequence of asymptotic cones Cg k (resp. C^ ar ' k ) 
converges to C% (resp. C£ ar )• We will use the following terminology [T7] : 

• The fatness ©(P^ -1 ) of a polyhedron p^ -1 is defined as follows: If a 
is an /-simplex, we begin to define the size e(p) of a: It is the maximum 
over all edges e of a of the length of e. Then, the fatness of a is the real 
number 

0 , x . . VOlj(fJb) 

B(cr) = mm mm — . 

fi ^-simplex in a j'e{l,...,Z} £{0~) J 

Finally, the fatness of P is the minimum of the fatness of its simplices. We 
denote by Fq the class of polyhedra in K N with fatness greater or equal 
to 0. 

• An (TV — 1) dimensional polyhedron p N ~ x in K N is closely inscribed in 
a smooth hypersurface W^ -1 if its vertices belong to W N ~ 1 and if the 
orthogonal projection of P N_1 onto W^ -1 is a bijection. 

• Let P N_1 be an (TV — 1) dimensional polyhedron (bounding a domain V) 
closely inscribed in a smooth hypersurface W N ~ 1 in E^. The angular 
deviation a m between m G P N ~ X and pr W N-i(m) is the maximal angle 
between the normal £ P r wN _ 1(rn) of W N ~ 1 and (m,n m ), where (m, n rn ) 

belongs to the support of N(X>). If B is any Borel subset of E^, we write 

olb = sup a rn . 

One of the classical observations is that a "good" fatness of a polyhedron 
and a "small" Hausdorff distance of this polyhedron to a smooth hypersurface 
in which it is closely inscribed imply that the angular deviation is "small" . 

4.1 A convergence result 

Let us state now our convergence theorem. For any cone C, we denote by UC 
the basis of C; that is, UC = Cn§ Ar_1 (0, 1), where S 7V_1 (0, 1) is the unit sphere 
centered at the origin. We also use the distance d defined on the class of subsets 
of E^ by 

d(A,B)= inf d(x,y). (36) 

xeA,yeB 

Theorem 1 Let D be a (compact) domain ofK N bounded by a smooth hyper- 
surface W N ~ 1 . Let (T>k) be a sequence of domains ofK N bounded by polyhedra 
(P^ -1 ) closely inscribed in W N ~ 1 such that: 
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1. The limit of (Pj^ 1 ) is W N 1 for the Hausdorff distance. 

2. The fatness of (P^ -1 ) is uniformly bounded from below by a non-negative 
constant: there exists 0 > 0 such that for all fcGN ; ^ V_1 £ 

Let G Sjgiv, siAcft /or all x G E^, |<£f)|(dl?) = 0. Then, every sequence 

(xk G C^ fe ) of unit vectors admits a subsequence (still denoted by (xk G C^ k )) 
that converges to a unit vector of C% . In particular, 

lim d(UC^ v \UC^ D ) = 0. (37) 

In the smooth case, the assumption on the boundary of B can be translated 
in terms of the second fundamental form of W N ~ X . We say that the normal 
curvature at a point p of a submanifold V of W N ~ X is null if the second funda- 
mental form h p of W N ~ X satisfies h p (u,u) = 0 for every u tangent to V at p. 
We get: 

Corollary 1 Under assumptions [7] and\^ of Theorem^ suppose that B is a 
Borel set such that W^ -1 D dB is smooth and with null normal curvature. 
Then, the conclusion of Theorem^ holds. 

To prove this corollary, we simply remark that under these assumptions, 

§ x D (dB) = / h(pr TW N-ix,pr TW N-ix) = 0. (38) 

Jw N ~ 1 r\dB 

The rest of this section focuses on the proof of Theorem [T] 
4.1.1 Convergence of sequences of normal cycles 

We need to introduce the flat norm on V t (E N ) as follows. The mass of an 
/-current T is the real number 

M(T) = sup{2»}, (39) 

where the supremum is taken over all /-differential forms with compact support 
such that sup^]]^ \ lj0 m\ < 1- The flat norm of an /-current T is the real number 

F(T) = inf{M(A) + M(B)}, (40) 

where the infimum is taken over all rectifiable /-currents A and (/ + l)-currents 
B such that T = A + dB. Our main ingredient in our study of convergence and 
approximation of the asymptotic cones is the following result, which is a simple 
reformulation of Theorem 67 of [17] (page 200) for polyhedra: 

Theorem 2 If p N ~ x is a closed (N — 1) dimensional polyhedron bounding a 
domain D and closely inscribed in a smooth closed hypersurface W N ~ 1 ofK N 
bounding a domain V, then for any Borel subset B of P N ~ X , 

<K(5 B +a B )MCN(D) lTBEN ), 

where 5b is the Hausdorff distance between B and pr W N-i(B) and K is a con- 
stant depending on the norm of the second fundamental form ofW N ~ 1 . 
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The following proposition can be deduced from Theorem [2] in a slightly 
different version, see also [T3] . 

Proposition 3 Under the assumptions [7] and^ of Theorem^ 

1. The masses M((N(2\))) are uniformly bounded from above; 

2. The sequence (N(P/ e )) converges to N(D) for the flat norm. 

4.1.2 Convergence of sequences of asymptotic measures 

Our framework is the space of (signed) Radon measures on K N with finite total 
variation, endowed with the norm ||.||i defined for every fi (with finite total 
variation |/i|) by 

IHIi = / EJV <M 

It is well known that this space is the (topological) dual to the space C C (E N ) of 
continuous functions with compact support on E^, endowed with the norm |.| 
defined by 

11/11= sup \f(x)\. 

xeE N 

The space of (signed) Radon measures (with finite total variation) can also be 
endowed with the topology of weak convergence of measures: a sequence of 
Radon measures (jik) on E N (weakly) converges to fi if, for every continuous 
function / with compact support on E N (resp. TE N ), jJLk(f) converges to fi(f). 

Proposition 4 Under the assumptions^ and\^ of Theorem^ 

1. For each vector x G E N , the sequence of measures ($f> ) converges to <& X D 
for the weak convergence of measures on K N ; 

2. For each unit vector x G K N , the sequence of measures (3>f> fc ) is ||.||i- 
bounded. 

Proof of Proposition [4} Item[T]is a direct consequence of Proposition [3J since 
the flat convergence of the sequence of normal cycles (N(Dfc)) implies the weak 
convergence of the measures (3>f> fc ). Item [2] is an application of the Theorem 
of Banach-Steinhaus: Since ) converges to § X D for the weak convergence 
of measures on E N , for each / G C C (E N ), sup fc | < f >< +00. The 

Theorem of Banach-Steinhaus then implies that ||3>f>J |i is uniformly bounded 
with respect to k; that is, the sequence of measures (|$x>J) ^ s ||-||i-bounded. 

Let us now explain our assumption on the boundary of the Borel subset B 
in Theorem [TJ In general, the weak convergence of the sequence (3>f> fc ) does not 
imply the convergence of (B)) to (^^(B)) for every Borel subset. Indeed, 
generically, characteristic functions are not continuous. That is why we restrict 
our study to a class of Borel subsets with suitable boundaries with respect to 
D, such that we can use the following general lemma (see [16 Chapter 4 for 
instance): 

Lemma 1 Let (fik) be a sequence of (signed) Radon measures on K N such that 
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1. (/ifc) converges to \i for the weak topology, 

2. the sequence (|/ifc|) of total variation of (fi k ) converges to a Radon measure 
v, for the weak topology. 

If the boundary dB of B G Sj^at satisfies v(dB) = 0 ; £ften 

lim fi k (B)=fi(B). (42) 

/c— ^oo 

Since the sequence (|3>f>J) * s ||-||i-bounded, we can extract a subsequence 
that converges. From Proposition [3] and Lemma [I] we deduce: 

Proposition 5 Under assumptions [7] and ^ of Theorem^ if x G E(E N ) and 
B G B e n satisfy \$ x D \(dB) = 0, then 

lim $% k (B) = ** D (B). (43) 

The last step of the proof of Theorem [T] is to relate the behavior of the 
sequence (3>f> ) f° r an y x t° tne behavior of their associated quadratic cones. 
For a fixed B G B-^n , we will study the quadratic forms x —> $f, fe (i?) and 
x — » ^^(i?) introduced in Section 3.2 We use the norm of uniform convergence 
on the space of quadratic forms: A sequence of quadratic forms (q k ) defined on 
K N converges to a quadratic form q if sup^n^ \q k (x) — q(x)\ tends to 0 when 
k tends to infinity. 

Lemma 2 Let (q k ) be a sequence of quadratic forms defined on K N , which 
converges to a quadratic form q. Let (x k ) be a sequence of unit vectors in 
E N , such that for each k G N, x k belongs to the isotropic cone C k of qk (i.e. 
Qk(%k) = 0,). Then there exists a subsequence of (x k ) that converges to a unit 
vector x belonging to the isotropic cone C of q. 

Proof of Lemma [2j We have 

lim sup \q(z) - q k {z)\ =0. (44) 

fc^oo| N | = 1 

Suppose that (x k ) is a sequence of unit vectors such that for all /c, x k G C k . 
Then, by the compacity of the unit sphere, there exists a subsequence (that we 



still denote by that converges to a unit vector x. From 44 , we deduce 



lim \q(x k ) - q k (x k )\ lim \q(x k )\ \q(x)\ 0, 

k^oo k^oo 

which means that x G C. 

We remark that under the assumption of Lemma |2j we cannot claim that 
the sequence of Hausdorff distances between UC k and UC tends to 0 when k 
tends to infinity, as shown in the following example. Consider the quadratic 
forms in E 3 defined for every k G N* by 

q k (u,v,w) = ^(u 2 + v 2 + w 2 ). (45) 

We have for every k G N*C fe = {0}, C = E N , UC k = 0 and UC = S N_1 (0, 1). 

The proof of Theorem [I] follows from Proposition [3J Lemma [I] Proposition 
[Hand Lemma [21 
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4.2 An approximation result 

In this section, for simplicity, we restrict our study to surfaces in E 3 . We assume 
that W 2 (bounding D) and P (bounding V) are fixed, with P being closely 
inscribed in W 2 . We suppose (without any restriction) that P is endowed with 
a triangulation, denoting by t a generic triangle, and by r(t) its circum-radius. 
The following theorem compares the asymptotic cone of V over a Borel set B 
composed of a union of triangles of P and the asymptotic cone of D over the 
orthogonal projection pr W 2(B) of B on W 2 . If C is the isotropic cone of a 
quadratic form q, we denote by A e (UC) the set of unit vectors x e-close to C; 
that is, \q(x)\ < e. 

Theorem 3 Let D be a domain o/E 3 bounded by a smooth hypersurface W 2 . 
Let V be a domain o/E 3 bounded by a polyhedron P G Tq, 6 > 0 closely inscribed 
in W 2 . For any e > 0, there exists r] > 0 such that if max{r(t),t G B} < n, 
then 

UC p**,v cA e (UC^£ (B) ). 

In other words, roughly speaking, under the assumptions of Theorem |3j the 
asymptotic cones of V are close to the asymptotic cones of D. The proof uses 
the results of Section 5-2 of [8]. We summarize them in the following proposition. 

Proposition 6 Under the assumptions of Theorem^ for every unit vector x, 

\$UB) ~ *f,(pr wa (B))| = Kmax{r(f),t e B}, (46) 

where K is a constant depending on the area of B, the length of its boundary, 
the geometry ofW 2 , and 0. 



If x is a unit vector belonging to UC^ r,T> , then 46 implies that 

| | =Kmax{r(t),tG B}. 
Consequently, if the triangles of P have a sufficiently small circumradius, then 

\** D (pr w2 (B))\<e. 

The conclusion follows. 



5 Some experiments 

5.1 Construction of asymptotic directions of a triangula- 
tion 

To mimic the smooth situation, it may be interesting to deduce asymptotic 
directions and asymptotic lines from the asymptotic cones defined on a singular 
surface. This construction may be achieved if one has a natural plane field 
on this surface. This is the case for triangulated surfaces, where each triangle 
spans a plane. Then, to build asymptotic directions at a point of a triangulated 
surface P 2 , associated to a Borel set £?, one can proceed as follows: 
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• Consider a point m on P 2 (bounding V) and build a Borel set B around 
m; for instance, a ball centered at m with a suitable radius. 

• Build the asymptotic cone C^ r,T> whose vertex is m, associated with B. 

• Build the intersection of C^ r,T> with the plane spanned by the face of 
the triangle that contains m. The result is two directions, called the 
asymptotic directions of P 2 at m. 

• When m runs over the surface, build a cross field (reduced to a point when 
the cone is reduced to a point). By integrating the cross field, we have 
asymptotic lines. 




(a) (b) 

Figure 4: (a) At each point where the asymptotic cone is not reduced to {0}, we 
build two asymptotic directions on the top of the Lilium tower (b) By integration 
we have asymptotic lines. 




(a) (b) 

Figure 5: (a) The cross field of asymptotic directions on the Heydar Aliyev 
Center in Baku (b) The asymptotic lines obtained by integrating the asymptotic 
directions 
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5.2 Approximation of the asymptotic lines of a smooth 
surface 



Using the construction of asymptotic cones, we can approximate the asymptotic 
directions (resp. lines) of a smooth surface W 2 in E 3 . We give an explicit 
example here. First of all, let us consider a portion of a (smooth) catenoid W 2 , 
and a triangulation P 2 closely inscribed in W 2 , with a sufficiently dense set of 
vertices. As shown in Figures 6a| 6bJ 6c , 6d , the intersection of the tangent plane 
of W 2 at a point m with the asymptotic cone of a ball centered at m is reduced 
to two lines that are a discrete approximation of the asymptotic directions of 
W 2 at m. If necessary, we can approximate the tangent plane itself by the plane 
spanned by a face of the triangulation. By integrating the directions field, we 
obtain discrete asymptotic lines, as shown in Figures 7a and 7b The reader 



can compare the asymptotic lines directly built on the smooth catenoid (these 
lines are orthogonal since the catenoid is a minimal surface) with the "discrete 
ones" obtained by this procedure. 





(c) (d) 

Figure 6: (a) Step 1: A green asymptotic cone built at a point m of a triangu- 
lated surface approximating a (smooth) catenoid (b) Step 2: The blue plane is 
an approximation of the tangent plane of the (smooth) catenoid at m (c) Step 3: 
The intersection of the green asymptotic cone and the blue plane gives two lines 
intersecting at m (d) Step 4: When m runs over the triangulation, we obtain a 
cross field that approximates the field of asymptotic directions of the (smooth) 
catenoid 
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(a) 



(b) 



Figure 7: (a) Integrating the cross field gives an approximation of the asymptotic 
lines of the catenoid (b) The cross field is an approximation of the field of 
asymptotic directions on the (smooth) catenoid. They can be compared with 
the blue lines, which are the asymptotic lines directly computed on the (smooth) 
catenoid 

5.3 Comparison of asymptotic lines 

We remark that the previous example gives (roughly speaking) "very good" 
results because the triangulated polyhedron approximating the catenoid suits it 
correctly, in the sense that the "good" fatness of the triangles implies that the 
normal of any triangle t is close to the normals of the orthogonal projection of t 
on the catenoid. To be more precise, we can estimate the error er as follows: Let 
W 2 be a smooth surface approximated by a triangulated polyhedron P 2 closely 
inscribed on it. Let m be a vertex of P 2 , and let us denote by e\ and e 2 the 
asymptotic directions of W 2 at m, and by ei, e 2 the approximated asymptotic 
directions at m. We define 

er m = inf (i(Z(ei, ei) + ^(e 2 , e 2 )), ^(Z(e l5 e 2 ) + Z(e 2 , ei))) , (47) 
where all the angles belong to (0, |). 

For instance, let us consider the portion of the Enneper surface shown in 
Figure [8j and a triangulated polyhedron closely inscribed on it. The error er is 
always less than or equal to 5 degrees. 
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(a) (b) 

Figure 8: (a) Comparison of discrete and smooth asymptotic lines on the En- 
neper surface (b) Some asymptotic lines computed on the polyhedron 

In the following example, we show that, in accordance with the theory, the 
error may be large even with a very thin triangulation closely inscribed in a 
smooth surface, if the angle between the tangent plane at a point and the 
corresponding triangle is too large (the same phenomenon appears when one 
compares the area of a cylinder with the area of a Lantern of Schwarz inscribed 
on it, see [17 for instance). Here, we consider a (smooth) surface W 2 in R 3 
obtained as the graph of a tensor product 5-spline function 

/ : G = [0, 1] x [0, 1] -»• R, 

of degree 3 in each variable, defined by 5 x 5 control points over a regular grid 
of 80 x 80 points. Each square of G is triangulated by taking both diagonals. 
We build the corresponding (piecewise linear) triangulation P 2 inscribed on 
W 2 . Then, we compare, at each vertex m of the triangulation, the normal of 
W 2 at m and the average of the normals of the triangles incident to m. The 



error varies between 0 and 0.5 degree (see Figure 9b). On the other hand, we 
compute on the same triangulation the asymptotic directions of W 2 and the 
asymptotic cones of P 2 , from which we deduce discrete asymptotic directions 
by intersecting the cones with the tangent planes as before. Then, we compare 
at each point m the error er m given in Figure [9a] The error at each vertex m 
varies between 0 and 5 degrees, according to the behavior of the normal of the 
faces incident to m with respect to the normal of W 2 at m. 
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(c) 

Figure 9: (a) Comparison of the normals computed on a 5-spline surface and 
an approximating triangulation (b) Comparison of the asymptotic directions by 
computing er m at each vertex m (the black points are convex points, where the 
asymptotic cone is reduced to {0} (c) Some asymptotic lines computed on the 
triangulation 



5.4 Deformation of asymptotic lines of discrete surfaces 

In the following example, we produce a deformation of "discrete" asymptotic 
lines as follows: We build a triangulation closely inscribed on a smooth surface 
W 2 (here a Chen's surface). We then compute asymptotic lines by the previous 
process using balls of radius R = 3 (the normalization is such that the average 
length of the edges is 1). Then, we slightly modify the position of the vertices 
that can now be out of W 2 (in other words, we create noisy data), without 
modifying the topology of the triangulation. With this new triangulation, we 
compute new asymptotic lines using the same process and using balls of the 
same radius (R = 3) or a different radius (R = 6). 



22 



(a) (b) (c) 

Figure 10: By perturbing the positions of the vertices of the approximated 
Chen's surface, we get noisy data and corresponding noisy discrete asymptotic 
lines (a) The initial triangulation with asymptotic lines (R = 3) (b) Asymptotic 
lines with noisy data and R = 3 (c) Asymptotic lines with noisy data and R = 6 

6 Conclusion and perspective 

In this article, with any smooth subspace or singular geometric subspace W of a 
Riemannian manifold M, we associate a family of cones, defined over any Borel 
subset of M . These cones are the generalization of the asymptotic directions 
defined at each point of a smooth surface of the Euclidean space. We obtain 
convergence and approximation theorems, when a sequence of polyhedra tends 
to a smooth subspace. As a consequence, we find good approximations of the 
asymptotic lines of a surface when it is approximated by a suitable triangulation. 
In our future work, we will study the relations between the characteristic of these 
cones and the geometry (and topology) of W. We will apply these results in 
different fields (face recognition to detect similarities between scanned faces for 
instance). 
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